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ABSTRACT 

Surface profiles were generated by a fractal algorithm and 
haptically rendered on a force feedback joystick. Subjects were 
asked to use the joystick to explore pairs of surfaces and report 
to the experimenter which of the surfaces they felt w ( as rougher. 
Surfaces were characterized by their root mean square (RMS) 
amplitude and th^ Ajcta! dimension. The most important 
factor affecting tnc perceived roughness ot the fractal surfaces 
was the RMS amplitude of the surface. When comparing 
surfaces of fractal dimension 1.2-1.35 it was found that the 
fractal dimension was negatively correlated with perceived 
roughness. 

INTRODUCTION 

The perception of texture is an important means by which 
humans identify surfaces around them. In field geology, 
scientists use texture to help identify rocks and determine their 
history (e.g., the amount of weathering they have been exposed 
to). The work in this paper is part of an effort to allow field 
geologists to explore remote planetary surfaces, displayed as 
virtual surfaces in an immersive haptic and visual display. An 
important component of the perceived texture of rock surfaces is 
their roughness, and the focus of this paper is the display of 
roughness using models that produce surface profiles similar to 
those found in nature. 

In particular, we have identified two parameters, RMS 
amplitude and fractal dimension, that are useful for generating 
synthetic surfaces similar to those measured using an optical 
profilometer on rock surfaces. In this paper we report on 
experiments to determine how human subjects relate perceived 
roughness to variations in these two parameters when 
interacting with the surfaces via a haptic display. 
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PREVIOUS WORK 

Earlier surface roughness perception experiments had 
subjects touch metal gratings consisting of rectangular grooves 
[10][ 12]. These studies showed that the dominant factor in 
determining perceived roughness was the groove width. More 
recent work by Lederman and Klatzky has argued for the 
viability of using a probe to encode vibration informat n .j 
discriminate roughness [8][9]. During these experiments 
subjects used a probe to explore surfaces made up of a pseudo- 
random pattern of raised dots. It was found that the perceived 
roughness of the patterns when using a probe or a finger 
increased with increased interelement spacing until reaching a 
peak and then decreased [8]. 

A common. 'feature of these experiments is that they use 
artificially constructed surfaces of predetermined heights. While 
the gratings are more deterministic than the spatially 
randomized raised dots, the dot patterns still have only two 
heights: the height of the dot and the flat surface it is resting on. 

Other work has examined the use of stochastic methods for 
texture display [16][ 17]. In the present case, we are interested in 
the perceived roughness of irregular virtual surfaces whose 
geometric properties (e.g., peak heights and spatial intensity) 
match those found on rocks. 

We selected a fractal technique to simulate the surface 
profiles. Previously we had found this technique to favorably 
mimic rock surface profiles when compared in terms of 
standard metrics from surface metrology, e.g. mean amplitude 
departure from a mean reference line. RMS amplitude deviation 
from a mean reference line, peak density, and kurtosis [5]. 



Characterizing Surface Profile Roughness with Fractal 
Dimension and RMS Amplitude 

Before synthesizing any profiles, the metrics used to 
characterize the surfaces must quickly be discussed. We use two 
parameters to describe the surfaces: the fractal dimension from 
Iractal geometry, and the RMS amplitude from traditional 
surface metrology. 

Fractals have been used to describe irregular shapes that do 
not lend themselves to description by Euclidean geometry. 
Natural structures such as mountains, coastlines, clouds and 
snowflakes, in addition to recursive, self-similar structures such 
as the Koch snow flake curve are examples of shapes that exhibit 
a fractal nature [13J. Each of these objects has a non-integer 
fractal dimension, measuring how much space they occupy. For 
example a straight line has a topological dimension of one, and 
a square has a dimension of two. The higher the dimension, the 
more space the curve occupies. Details of measuring the fractal 
dimension can be found in The Science of Fractal Images [13]. 
Several attempts have been made to use the fractal dimension to 
characterize the roughness of rock surface profiles [6][7]. A 
fractal analysis of a surface can yield a fractal dimension and an 
amplitude coefficient [ 13][3]. 

Amplitude measures have been a popular way in surface 
metrology to describe the surface roughness of a profile [5]. The 
amplitude measure that we use is R q . the root mean square 
deviation from the mean line. In continuous form R ( is defined 
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where L is the length of the profile. When measuring fractal 
properties of real surfaces we were able to compute the RMS 
deviation from the fractal amplitude coefficient. With this 
framework we can begin to synthesize rough surfaces. 

Synthesizing a Surface: Fractal Motion using Fourier Filter- 
ing 

For our experiment we desired to generate surfaces 
mimicking those found in nature, while being able to control 
their frequency and amplitude characteristics. Fractal geometry 
has been used in computer graphics to generate mountainous 
landscapes, plants, and water surfaces. Using an appropriate 
algorithm we can vary the fractal dimension and the root mean 
square amplitude of our generated surfaces. The first algorithm 
we tried was based on the "Fractal Motion using Fourier 
Filtering method" [I3J. Fourier series has been used before in a 
similar algorithm to generate surfaces [18]. The Fourier 
Filtering method differs by imposing the power spectral density 

condition of fractal Brownian motion, P(J)a(\/f , onto the 
coefficients of the discrete Fourier transform. However this 


condition uses only the fractal dimension. This does not allow 
adjusting of the amplitudes of the generated surfaces. 

Ganti uses a similar I //‘power spectral density function to 
generate fractal surfaces but also includes an amplitude 
parameter C that scales the amplitudes of frequencies [3| 
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where P is the power spectral density at /rth frequency/, D is the 
fractal dimension, C is the amplitude coefficient, L is the length 
of the profile, and r| is the instrument measurement resolution 
of the profile. For synthesizing profiles we used the encoder 
resolution for r\ . 

When measuring profiles, the amplitude coefficient C can be 
related to the RMS amplitude of the heights by 

V (jf ~- r o) = Cq- D " V'- D (2 ) 

where t is the sampling resolution [3]. For synthesizing profiles 
we took t and rj to be same. This equation enables us to 
synthesize a surface with any desired root mean square 
amplitude. 

To synthesize a surface the conditions of equation (1) are 
imposed onto the coefficients of the discrete Fourier transform 
[I3][3] 
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[31 tor a desired profile of A/ points, random sign is randomly ± 1 
and rand is a random Gaussian number from 0 to 1. We then 
form H\ using the complex conjugate operator * 
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The conditions of the power law are imposed onto H k by 
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where/is the fundamental frequency. 1/L, the inverse of the 
desired length of the profile. We then compute our synthesized 
function. X(t ), using our random coefficients and the discrete 
inverse Fourier transform [13] 


Xd) = V (Rc!: t ' (h!) + ImH k sin ( hi ) ) 
k = 1 

One problem with using this method is that the generated 



Figure 1. 100 Grit Painter’s Sandpaper Profile (top). 
Fractal Simulations D=1.166 and C=.0063 (middle), 
C=.12 (bottom). 


surface repeats itself every 2rtn intervals. To circumnavigate 
this problem we divide the surface we would like to generate 


into 33)552 mm fragments of 128 equally spaced points, 
.generating a new set of Fourier coelt.c.ent pairs lor eac 
fragment. The newly generated (ragment has a constant added 
to u so that it IS attached to end ot the last tragment. This 
process continues until the new surface is complete 

In practice we also found that when generating profiles with 
amplitude coefficients based on equation (2). the profiles had 
R s less than predicted. We increased the amplitude coefficient 
C until a surface was generated with the desired RMS 
amplitude. During our previous experiments when we ana yze 
actual surfaces, the measured amplitude coefficient and trac a 
dimension did predict the measured R q . 

A real profile of 100 grit painter s sandpaper and two fractal 
simulations based on its measured fractal parameters D and C. 
are shown in Fig. 1. The middle profile exhibits the algorithm s 
problem with producing surfaces with amplitudes smaller than 
that predicted by equation (2). The last profile had an increased 
amplitude coefficient C to increase its R q to approximately 
same R q as the painter's sandpaper. 

EXPERIMENTAL SETUP 

Fractallv simulated surface profiles were haptically rendered 
on an Immersion Impulse Engine 2000 force feedback joystick. 
Fourteen subjects, eleven male and 3 female, from ages - I to 
34 participated in ihe experiment. Subjects were naive to e 
purpose of the experiment All subjects had a- least minima! 
experience using haptic devices. SuDje^ts w. e -»ked to use t e 
joystick to explore pairs of surfaces and report to e 
experimenter which one they felt was rougher, surface one or 
surface two. One surface was displayed at a time. Subjects 
controlled which surface was displayed by clicking the joystick 
trigger. Subjects were allowed to switch between surfaces as 
often as they desired and to explore surfaces as long as they 
desired. Sixteen pairs of surfaces were given to subjects in 
random order. Assignment of surfaces as surface one or two was 
also randomized. Typically subjects spent fifteen to twenty 
minutes exploring surfaces, with the shortest time being 
approximately ten minutes and the longest being approximately 
thirty-five minutes. 

Two different reference surfaces were compared to eight 
surfaces each to form the sixteen trial pairs. Each reference 
surface was compared to surfaces with a combination of lesser, 
equal, and greater fractal dimension and RMS amplitude. Figure 
o demonstrates the surface group pairings. Reference surface 
one is compared against each of the surfaces of group one while 
reference surface two is compared against the group two 

surfaces. 

The surfaces used in the experiment were generated in 
lengths and patched together to form 80mm lengths. Eac 
length was offset to meet that last one so that no discontinuity 
between lengths could be detected. Surface profiles were 






resampled to a resolution of 0.!545mm/sampl c . five times the 

encoder resolution of the invalid- Th. . 

1 me joystick. I his saved memory without 
compromising the feel of the surface. 

aboveT P,mg K thC SU , rtaCe alS ° SCtS " S Spa " al rcso,ution tojust 
The he T " bandw,dth of 'he force feedback joystick 

of iCr?? f °T ° f UP 8 9 N Wi,h a bandwid 'h 
. f lth J iypita minimu m user lateral exploration 

■ peed of approximately 20 mm/s. the achievable^ spatial 
resolution considering the bandw.dth is about 0.17 cycles/mm. 

Contact Model 

Subjects explore the virtual surfaces through the force 

contoct Vilely US,n§ h dy " amiC m ° dd ° f S ‘ y,US t0 surface 
contact, the forces are d, splayed to a user resultmg from 

interaction w„h the surface profile data. Previously !- e had 

begun with an initial model similar to the Sandpaper svstem 

heiihts tangemia ' f0rce Redback based on the change in 
eights with a vertical restoring force based on the penetration 

to be uns!t f Pr ° ^ " Pr£ViOUS Pr ° ject we found 'his model 

rocks m wlTV 'T erS ' ng SUrf3Ce Pr0fl,es taken from 
rock, fl] . We therefore implemented a dynamic model (Fig 3) 

at represents the normal and tangential forces at the stylus- 

sur ac e contact and accounts for the possibi)jty 

the surface 0nt This and ^“ nCmg ° r skippin g over the valleys of 

collected bv W3S USCd t0 dlSp,3y r ° ck ^faces 

dunng a NASA T SenSOr t0 planetar > geoiog.sts 

... g ‘ ^ Ames field experiment in February 1999 m 

, r : be llr nt described in th.s paper «e used this mode! 
haptically display the synthesized fractal profiles. 

First a height profile as a function of position, v (x) is 
generated by our fractal algorithm. The stylus is then modeled 
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Figure 2. Surface Croup Parings 


as a mass, spring, damper system connected to the vertical mm.t 

° ‘ U “ rS ,P,,t m,ertace devicc - V The horizontal position 
of the stylus is dtrectly coupled to the horizontal p„s,„on of the 
interlace. While in contact w„h the surface the stvlu n 1 
are computed by " > n amics 

= A t v 'r. v /0) + rf| r v^.v ( /)) + f 

For our application the parameters were tuned for feel to eve 
our virtual stylus parameters of K=\ .28 N/mm and negl.g.ble M 

as?fncn n n a r ^ ° f ^ S,y ' US aga, " St the surface ls modeled 
a frictionless contact point. As illustrated by the magnified 

porno„ of F, g . 3. the no™, f„ rcc , s we „d,cnJ,o“ 

tangent of the contact point. It is this reaction force that ,s 
isp aye to the user. The tangent is computed by taking the 
do„v a „ve of Ihe surface with respec, ,o the x coordinate J,et 
y In “"-Wlicedon this derivative is precomputed, as we' 
know the surface profile a priori. 

The normal force is the sum of x and y component forces. 

F n = K * K 

The y component represents the vert, cal react, on force The x 
component represents the horizontal react, on force 


F r = /. 
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If the user's position y, is ever greater than the profile height v 
contact ,s broken with the surface and contact forces are set to ' 
zero. This allows the user to bounce off and back onto the 
surface, especially while moving fast. 

Because the subject can expk.-e the surface freelv. exact force 
evels displayed to the subject are dependent on how thev chose 

urfa P H "I™' C g - dCpth ° f Penetrat - ‘"to the v^nual 
surface, and speed of exploration. 

feil Jlass^PlThT ** * UrfaCCS modeled without friction 
reel glassy [-]. While we did not use friction in this model 

Str** r e of ,he subj “ ,s 

fe ‘ glaSSy ° r Sllpper y There is a small amount of 
friction inherent to the force feedback device we used The 

Impulse Engine 2000 specifications report a maximum 
backdnve friction of 0.14 N [14], Its possible that a 
combination of factors including the forces displayed due o the 
irregular, ty of the surface, the backdrive fnction and device 

friction mmateS a " y S ‘ iPPery fed by Creating an il,usion of 

RMS J !mn, V , Cl H mteraCti ° ns with surfaces with smaller 

RMS amplitudes ( 2mm-.27mm, usmg this model feels like 

usmg a stylus to stroke sandstone, while m.eractmg with fra-f.l 
surfaces w„h larger RMS amplitudes ,>.5mm) , i k 





stroking a weathered conglomerate rock or broken concrete 
surface. 


DATA AND ANALYSIS 

Figures 4 and 5 show the perceived roughness of a surface 
when compared to a reference surface. Data points are plotted 
along the x and y axis by their fractal dimension and RMS 
amplitude. The center spot labeled ref represents the surface that 
was compared to the other eight surfaces surrounding it. The 
percentage of subjects w'ho reported a surface as rougher than 
the center point surface is plotted along the z-axis. For example 
in Fig. 4, / j c ot subjects reported that the surface of fractal 
dimension 1.35, RMS amplitude 0.3502mm. was rougher than 



Figure 4. Perceived Roughness vs. Surface D=1.275, 
RMS Amplitude=.2781 



Figure 5. Perceived Roughness vs. Surface 
D=1.425, RMS Amplitude=.75 

the reference surface of dimension 1.275 and RMS amplitude 
0.278 1mm. Likewise the other data points in the bar graph are 
of the percentage perceived rougher when compared to the 
center surface of fractal dimension 1.275, RM,, amplitude 
0.2781mm. The reference surface response is set to 50%. This 
would represent the case ot comparing the reference surface to 
itself, although for fatigue considerations this case was not 
actually presented to subjects. 

Looking at Fig. 4 and 5 across the RMS amplitude axis, both 
cases show that surfaces with a higher RMS amplitude are 
perceived as rougher. For all cases of different RMS amplitudes 
between 71.4% to 100% of subjects reported the surface with 
the higher amplitude to be rougher. The mean percentage of 
subjects reporting that a surface with a higher amplitude to be 
rougher was 91.7% with a standard deviation of ±8.5% for both 
cases. 

It is much more difficult to discern any trends in roughness 
perception due to fractal dimension by only examining the data 
plots. To check for the statistical significance of the dimension 
and amplitude parameters we performed a logistic regression 
analysis. A logistical regression analysis was chosen because 
subject responses are binary. After a subject indicates which 
surface they believe is more rough, the result is recorded as 
positive or negative depending if they chose the test surface 
(positive) or the reference surface (negative!. 

For each case we fitted the subject responses to a multiple 
logistic regression model of the form [4] 
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using a software statistical analysis package [11]- Each 
coefficient (f, of the model represents a different experimental 

variable: fractal dimension, RMS amplitude, or a subject 
variable. After computing the coefficients and their respective 
estimated standard errors, we performed the Wald test for 
logistic regression models to check for statistical significance 
[4f The Wald test is conducted by first computing the univariate 
Wald statistic W = 3/(S£(P,)). where 5£(P,) is the standard 

error of the rth coefficient. Next the two tailed p-value of the 
Wald statistic is computed from the Chi-square on one degree of 
freedom distribution with a significance of a =.05. For any p- 
value less than .05 we consider the variable to be statistically 
significant. 

The p-values for the coefficients are presented in Table 1 and 
Table 2 on page 6. As expected for both cases the p-value of the 
Table 1: P-Valuesfor Perceived Roughness when Comparing to 
Reference Surface 1: D—l.275, RMS Amplitude—. 2 781 mm 


Coefficient Name 

Wald Statistic 

p-value 

Fractal Dimension 

-2.48803 

.0128 

RMS Amplitude 

4 47037 

7.8085e-6 

Subject 1 

-1.027 

.3044 

Subject 2 

- 1 .027 

.3044 

Subject 3 

1.027 

.3044 

Subject 4 

- 1 .5e- 16 

1 

Subject 5 

-.051 

.9593 

1 Subject 6 

-1.815 

.0696 

Subject 7 

-1.027 

.3044 

Subject 8 

-1.027 

.3044 

Subject 9 

-1.815 ' ~ 

.0696 

Subject 10 

-1.027 

.3044 

Subject 1 1 

-1.815 

0696 

Subject 12 

1 .027 

.3044 

Subject 13 

-1.027 

.3044 


RMS amplitude is well below 0.05, indicating that this variable 
is statistically significant. As the Wald statistic is positive in 
both cases, the likelihood of perceiving one surface rougher 
than the reference surface increases with increasing amplitude. 
For the second case of higher fractal dimensions. 1.35 to 1.5, the 
p-value for fractal dimension is 0.446. This indicates that this 


variable was not a significant (actor in perceiving surface 
Table 2: P- Values for Perceived Roughness when Comparing to 
Reference Surface 2: D=l J25, RMS Amplitude =.7 5mm 


Coefficient Name 

Wald Statistic 

p-value 

Fractal Dimension 

-.76212 

446 

RMS Amplitude 

4.9217 

8.58 le-7 

Subject l 

-6.63c-6 

l 

Subject 2 

-1.734 

0829 

Subject 3 

-1.734 

.0829 

Subject 4 

-2.504 

0123 

Subject 5 

-.966 

.334 

Subject 6 

-.966 

.334 

Subject 7 

-.966 

.334 

Subject 8 

-.966 

.334 

Subject 9 

-.966 

.334 

Subject 10 

-.966 

.334 

Subject 1 1 

-6.63e-6 

1 

Subject 12 

-.966 

.334 

Subject 13 

-719 

4721 


roughness when comparing surfaces of fractal dimension 
greater than 1.35. Interestingly in the first case when comparing 
surfaces of lower fractal dimension. 1.2 to 1.35, the fractal 
dimension parameter has a p-value of .0128 and a negative 
•Vaid statistic. Witn this p-value this vanaoie does meet our test 
for statistical signifance. This suggests that when comparing 
fractal surfaces with dimensions between 1.2 and 1.35, the 
surface with the lower fractal dimension contributes to it being 
perceived as rougher. 

This dichotomy between the higher and lower fractal 
dimension cases might exist because az the fractal dimension 
decreases the surface becomes more coarse, enforcing the 
perception of roughness. The difference between two surfaces 
with higher fractal dimensions may not be noticeable to a 
person. This is supported by Fig. 5. Examining the two trials 
where the RMS amplitude is the same, the percentage of 
subjects perceiving the surface of dimension 1.5 and the surface 
of dimension 1.35 as more rough than the reference surface of 
dimension 1.425, is 50% and 57%, respectively. This indicates a 
random preference. Subjects reported that these two cases were 
difficult to distinguish. 

CONCLUSIONS 

From examining both the data plots and the p-values the most 
important factor effecting the perceived roughness of fractal 
surfaces is the RMS amplitude. The logistic regression indicated 
that when comparing fractal surfaces with dimension 1.2-1.35, 
lower fractal dimension contributes to the perception of 
roughness. This agrees with previous work that reported 
roughness perception increased with increasing groove width 



and increasing interclement spacing [ tO|[ ! 2 j| SJ. The changes in 
those surface parameters would have decreased the fractal 
dimension as well. 

This effect does not hold w hen comparing surfaces of higher 
fractal dimensions - those over 1.35. For these comparisons 
fractal dimension was statistically insignificant. When 
comparing surfaces of the same RMS amplitude in the set of 
surfaces with fractal dimension 1.35 to 1.5, the subjects' 
selection were apparently random. 

In any case of fractal dimension. RMS amplitude is the 
overriding factor in determining surface roughness perception. 
When comparing surfaces of different RMS amplitudes, 
subjects selected the surface with the higher amplitude as 
rougher 71.4% to 100% of the time. 

ACKNOWLEDGMENTS 

The authors would like to thank Philip Beineke from the 
Statistical Consulting staff at Stanford University for his helpful 
advice on logistical regression and applying it to our data. 

REFERENCES 

[ 1 JOkamura, Costa, Turner, Richard, Cutkosky. Haptic Surface 
Exploration, ISER 99 

[2JGreen. D.F., ’’Texture Sensing and Simulation Using the 
PHANToM: Towards Remote Sensing of Soil Properties." in 
The Second PH ANToM User’s group workshop. MA. Oct 19- 
2 ~, 7 . 

[3] Ganti, S., Bhushan. B., Generalized fractal analysis and its 
applications to engineering surfaces. Wear, Vol. 180 pp 17-34 
1995. 

[4] Hosmer, D. W Lemeshow\ S.. Applied Logistic Regression. 
John Wiley & Sons. 1989. 

[5] Whitehouse, D. J.. Handbook Of Surface Metrology. Bristol ; 
Institute of Physics Pub.. 1994. 

[6] Jones, J.G.. Thomas, R. W., Earwicker, P. G., Fractal 
Properties of Computer Generated and Natural Geophysical 
Data. Computers and Geosciences. Vol. 15, No. 2. pp 227-235 
1989. 

[7] Brow n, S., A Note on the Description of Surface Roughness 
Using Fractal Dimension. Geophysical Research Letters, Vol. 

14. No. 11, pp. 1095-1098. November 1987. 

[8] Lederman. Klatzky, Feeling Through a Probe. DSC- Vol. 64, 
Proceedings of the ASME Dynamic Systems and Control 
Division, ASME 1998. 

[9] K!atzky, Lederman. Roughness Via a Rigid Probe: Effects of 
Exploration Speed. DSC- Vol. 67, Proceedings of the ASME 
Dynamic Systems and Control Division, ASME 1999. 

[ lOJLederman, S.J.. Tactile roughness of grooved surfaces: The 
touching process and effects of macro- and micro surface 


structure. Perception and Psychophysics. 32, pp. 109-116. 1974 

( 1 1 [Smyth, G. K. ( 1999). StatBox 4 0: A statistics toolbox for 
Matlab. University of Queensland, Brisbane 
(www. maths.uq.edu.au/-gks/matlab/statbox.html). 

( 1 2 [Taylor, M.M., & Lederman. S.J. Tactile roughness of 
grooved surfaces: A model and the effect of friction. Perception 
and Psychophysics, 17, pp. 23-36 1975. 

[ 1 3[Peitgen, H.-O., Saupe, D., The Science of Fractal Images. 
Springer- Verlah, New York, NY. 1988. 

[14]Impu!se Engine 2000 Specifications. http://w w w.force- 
feedback.com/research/research.html 

[15JM. Minsky, M. Ouh- Young, O. Steele. F.P. Brooks jr., M. 
Behensky, Feeling and Seeing: Issues in Force Display, Proc. 
Symp. Interactive 3D Graphics, pp. 235-243, 1990. 

[I6JJ.P. Fritz, K.E. Bamer, Stochastic Models for Haptic 
Texture, Proc. SPIE Ini. Symp. Intelligent Systems and 
Advanced Manufacturing Vol. 3901, pp. 34-44, 1996. 

[17] J.Siira, D.K. Haptic Rendering - A Stochastic Approach, 
Proc. 1996 IEEE Inti. Conf. Robotics and Automation, 
Minneapolis, Minnesota, pp. 557-562, 1996. 

[18] S.A. Wall, W.S. Harwin. Modeling of Surface Identifying 
Characteristics using Fourier Series. DSC- Vol. 67, Proceedings 
of the ASME Dynamic Systems and Control Division, pp. 65- 
71. ASME 1999. 




.*?< 0 = 0 + P| »| + + ■ + fiptp 

,V(l) 

\t 1 I 1 V 


using a software statistical analysis package [11]. Each 
coefficient of the model represents a different experimental 

variable: fractal dimension, RMS amplitude, or a subject 
variable. After computing the coefficients and their respective 
estimated standard errors, we performed the Wald test for 
logistic regression models to check for statistical significance 
[4]. The Wald test is conducted by first computing the univariate 
Wald statistic W = p./(S£((3 f )) , where SE($ t ) is the standard 

error of the ith coefficient. Next the two tailed p-value of the 
Wald statistic is computed from the Chi-square on one degree of 
freedom distribution with a significance of a =.05. For any p- 
value less than .05 we consider the variable to be statistically 
significant. 

The p- values for the coefficients are presented in Table 1 and 
Table 2 on page 6. As expected for both cases the p-value of the 
Table 1: P-Values for Perceived Roughness when Comparing to 
Reference Surface 1: D-1.275, RMS Amplitude ^27 8 1mm 


Coefficient Name 

Wald Statistic 

p-value 

Fractal Dimension 

-2.48803 

.0128 

RMS Amplitude 

4 47037 

7.8085e-6 

Subject 1 

-1.027 

.3044 

Subject 2 

-1.027 

.3044 

Subject 3 

1.027 

.3044 

Subject 4 

- 1 5e- 16 

1 

Subject 5 

-.051 

.9593 

i Subject 6 

-1.815 

0696 

Subject 7 

-1.027 

.3044 

Subject 8 

-1.027 

.3044 

Subject 9 

-1.815 ' ~ 

.0696 

Subject 10 

-1.027 

.3044 

Subject 1 1 

-1.815 

.0696 

Subject 12 

1.027 

.3044 

Subject 13 

-1.027 

.3044 


RMS amplitude is well below 0.05, indicating that this variable 
is statistically significant. As the Wald statistic is positive in 
both cases, the likelihood of perceiving one surface rougher 
than the reference surface increases with increasing amplitude. 
For the second case of higher fractal dimensions, 1.35 to 1.5, the 
p-value for fractal dimension is 0.446. This indicates that this 


variable was not a significant factor in perceiving surface 
Table 2: P- Values for Perceived Roughness when Comparing to 
Reference Surface 2: D- 1.425 , RMS Amplitude^ 5mm 


Coefficient Name 

Wald Statistic 

p- value 

Fractal Dimension 

-.76212 

446 

RMS Amplitude 

4.9217 

8.58 le-7 

Subject 1 

-6.63c-6 

1 

Subject 2 

-1.734 

.0829 

Subject 3 

-1.734 

.0829 

Subject 4 

-2.504 

.0123 

Subject 5 

-.966 

334 

Subject 6 

-.966 

.334 

Subject 7 

-.966 

.334 

Subject 8 

-.966 

.334 

Subject 9 

-.966 

.334 

Subject 10 

-.966 

.334 

Subject 1 1 

-6.63e-6 

1 

Subject 12 

-.966 

.334 

Subject 13 

.719 

4721 


roughness when comparing surfaces of fractal dimension 
greater than 1.35. Interestingly in the first case when comparing 
surfaces of lower fractal dimension, 1.2 to 1.35, the fractal 
dimension parameter has a p-value of .0128 and a negative 
'.Vaid statistic. With this p-value this vanaoie does meet our test 
for statistical signifance. This suggests that when comparing 
fractal surfaces with dimensions between 1.2 and 1.35, the 
surface with the lower fractal dimension contributes to it being 
perceived as rougher. 

This dichotomy between the higher and lower fractal 
dimension cases might exist because aJ the fractal dimension 
decreases the surface becomes more coarse, enforcing the 
perception of roughness. The difference between two surfaces 
with higher fractal dimensions may not be noticeable to a 
person. This is supported by Fig. 5. Examining the two trials 
where the RMS amplitude is the same, the percentage of 
subjects perceiving the surface of dimension 1.5 and the surface 
of dimension 1.35 as more rough than the reference surface of 
dimension 1.425, is 50% and 57%, respectively. This indicates a 
random preference. Subjects reported that these two cases were 
difficult to distinguish. 

CONCLUSIONS 

From examining both the data plots and the p- values the most 
important factor effecting the perceived roughness of fractal 
surfaces is the RMS amplitude. The logistic regression indicated 
that when comparing fractal surfaces with dimension 1.2-1.35, 
lower fractal dimension contributes to the perception of 
roughness. This agrees with previous work that reported 
roughness perception increased with increasing groove width 





unci increasing inlerclcinent spacing 1 10)1 12||8|. The changes in 
Ihnse surface paiamelers would have decreased Ihe Iraclal 
dimension as well. 

This effect does nor hold when comparing surfaces of higher 
Iraclal dimensions - ihose over 1.35. For fhese comparisons 
fractal dimension was statistically insigm ican . 
comparing surfaces of the same RMS amplitude in ‘he set of 
surfaces with fractal dimension 135 to 15, the subjec 
selection were apparently random. 

In any case of fractal dimension. RMS amplitude is the 
overriding factor in determining surface roughness perception. 
When comparing surfaces of different RMS amplitudes, 
subjects selected the surface with the higher amplitude as 
rougher 71.4% to 100% of the time. 
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